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Abstract 

We prove an equality-in-law relating the maximum of GUE Dyson's Brownian 
motion and the non-colliding systems with a wall. This generalizes the well known 
relation between the maximum of a Brownian motion and a reflected Brownian mo- 
tion. 



1 Introduction and Results 

Dyson's Brownian motion model of GUE (Gaussian unitary ensemble) is a stochastic pro- 
cess of positions of m particles, X(t) = (Xi(t), . . . , X m (t)) described by the stochastic 
differential equation, 

dt 

dXi = dBi+ V — , l<z<m, (1.1) 

i T~7 X ~ -^3 

l<j<m J 

where B i: l < i < m are independent one dimensional Brownian motions [5]. The process 
satisfies X\(t) < X 2 {t) < ■ ■ ■ < X m (t) for all t > 0. We remark that the process X can be 
started from the origin, i.e., one can take X((0) = 0, 1 < i < m. See [8]. 

One can introduce similar non-colliding system of m particles with a wall at the origin [6, 
7, 14]. The dynamics of the positions of the m particles X^ = (x[ , . . . , Xm) satisfying 
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< Xi(t) < X 2 (t) < ■ • • < X m (t) for all t > are described by the stochastic differential 
equation, 

dX i G) = dB i + ^5)+ E ( Y (0 1 Y (0 + Y (0 | Y (o )^ (1-2) 

X i l<j<m V X i - X j X i + X j / 

This process is referred to as Dyson's Brownian motion of type C . It can be interpreted as 
a system of m Brownian particles conditioned to never collide with each other or the wall. 

One can also consider the case where the wall above is replaced by a reflecting wall 
[7]. The dynamics of the positions of the m particles X^ = (x[ D \ . . . ,xffl) satisfying 
< Xi(t) < X 2 (t) < • • • < X m (t) for all t > 0, is described by the stochastic differential 
equation, 



dX^ 



dB i + U {i=1) dL(t)+ ( V (D) 1 V (D) + V (D) - V (D) )dt,l<i<m, (1.3) 

l<j<m \ X i ~ X j X i + X j J 



where L(t) denotes the local time of x[ D ^ at the origin. This process will be referred to 
as Dyson's Brownian motion of type D. Some authors consider a process defined by the 
s.d.e.s (1.3) without the local time term. In this case the first component of the process is 
not constrained to remain non-negative, and the process takes values in the Weyl chamber 
of type D, {\xi\ < x 2 < x 3 . . . < x m }. The process we consider with a reflecting wall is 
obtained from this by replacing the first component with its absolute value, with the local 
time term appearing as a consequence of Tanaka's formula. 

It is known the processes can be obtained using the Doob ^-transform, see [6]. 

Let (Pt' (C ' D) ;t > 0) be the transition semigroup for m independent Brownian motions 
killed on exiting {0 < X\ < x 2 ■ ■ ■ < x m }, resp. the transition semigroup for m independent 
Brownian motions reflected at the origin killed on exiting {0 < X\ < x 2 . . . < x m }. From 
the Karlin-McGregor formula, the corresponding densities can be written as 

det{0 t (xj - x'j) - (f) t (xi + 4)}i<*j<m, (1-4) 

resp., 

det{(j) t (Xi - Xj) + 4> t (Xi + 4)}l<iJ<m, ( L5 ) 

where <j) t {z) = -^e~ z ^ 2t \ Let 

m 

h^(x)=i[x, n {x]- x *i 

i=l l<i<j<m ^ Q'j 

hW( X )= n (^-^)- 

l<i<j<m 
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For notational simplicity we suppress the index C, D for the semigroups and in h in the 
following. Then one can show that h(x) is invariant for the P t ° semigroup and we may 
define a Markov semigroup by 

P t (x, dx) = h{x')P?{x, dx')/h(x). (1.7) 

This is the semigroup of the Dyson non-colliding system of Brownian motions of type C 
and D. Similarly to the X process, the processes and X™ can also be started from 
the origin (see [9] or use Lemma 4 in [7] and apply the same arguments as in [8]). 

In GUE Dyson's Brownian motion of n particles, let us take the initial conditions to be 
Xj(0) = 0, 1 < i < n. The quantity we are interested in is the maximum of the position of 
the top particle for a finite duration of time, max < s < i X n (s). In the sequel we write sup 
instead of max to conform with common usage in the literature. Let m be the integer such 
that n = 2m when n is even and n = 2m — 1 when n is odd. Consider the non-colliding 
systems of X^ C \X^ of m particles starting from the origin, XI C ' D) (0) =0,1 <i < m. 

Our main result of this note is 

Theorem 1. Let X and X^ C \X^ start from the origin. Then for each fixed t > 0, one 
has 

4 fx£\t), forn = 2m, 
sup X n (s) = t {D) (1.8) 
o<s<t I Am '(t), forn = 2m—l. 

To prove the theorem we introduce two more processes Zj and Yj. In the Z process, 
Z\ < Z 2 < ■ ■ ■ < Z n , Zi is a Brownian motion and Zj + i is reflected by Zj, 1 < j < n — 1. 
Here the reflection means the Skorokhod construction to push Zj + % up from Zj. More 
precisely, 

Z 1 (t) = B 1 (t), 

Z 3 {i) = sup (Zj^( S ) + Bj(t) - Bj(s)), 2 < j < n, ( L9 ) 

0<s<t 

where Bi,l < i < n are independent Brownian motions, each starting from 0. The process 
is the same as the process [X\(t), A~f (£), . . . , X£(t); t > 0) studied in section 4 of [15]. The 
representation (jl.9p was given earlier in [2]. In the Y process, < Yi < Y 2 < . . . < Y n , 
the interactions among 1^'s are the same as in the Z process, i.e., Yj + i is reflected by Yj, 
1 < j < n — 1 , but Y\ is now a Brownian motion reflected at the origin (again by Skorokhod 
construction). Similarly to (11. 9ft . 

Y 1 (t) = B 1 (t)- inf B 1 (s)= sup (B 1 (t) - B^s)), 

0<s<t 0<s<f 

- - (1.10) 

Yj(t) = sup (Yj^(s) + Bj{t) - Bj{s)), 2 < j < n. 

0<s<t 

From the results in [4,8, 15], we know 

(X n (t);t>0) = (Z n (t);t>0) (1.11) 
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and hence 

sup X n (s) = sup Z n (s). (1.12) 

0<s<t 0<s<t 

In this note we show 
Proposition 2. The following equalities in law hold between processes: 

(Y 2m (t);t>0)^(X^(t);t>0), 

(Y 2m _ 1 (t) ] t>o) = (x£\ty,t>o), 

m G N. 

The proof of this proposition is given in Section [21 The idea behind it is that the processes 
(Fj)i>i and (Xj D )j>i could be realized on a common probability space consisting of 
Brownian motions satisfying certain interlacing conditions with a boundary [15,16]. Such a 
system is expected to appear as a scaling limit of the discrete processes considered in [3,16]. 
In this enlarged process, the processes Y n (t) and Xm' D \t) just represent two different ways 
of looking at the evolution of a specific particle and so the statement of Proposition 2 follows 
immediately. Justification of such an approach is however quite involved, and we prefer to 
give a simple independent proof. See also [4] for another representation of Xm in terms 
of independent Brownian motions. 

Then to prove (11.81) it is enough to show 

Proposition 3. For each fixed t we have 

sup Z n (s)=Y n (t). (1.14) 

0<s<t 

This is shown in Section [31 For n = 1 case, this is well known from the Skorokhod 
construction of reflected Brownian motion [10]. The n > 1 case can also be understood 
graphically by reversing time direction and the order of particles. This relation could also 
be established as a limiting case of the last passage percolation. In fact the identities in 
our theorem was first anticipated from the consideration of a diffusion scaling limit of the 
totally asymmetric exclusion process with 2 speeds [1]. 
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2 Proof of proposition [2] 



In this section we prove the relation between X^ C ' D ^ and Y, (11.131) . The following Lemma 
is a generalization of the Rogers-Pitman criterion [11] for a function of a Markov process 
to be Markovian. 

Lemma 4. Suppose that {X(t) : t > 0} is a Markov process with state space E, evolving 
according to a transition semigroup (Pf, t > 0) and with initial distribution fi. Suppose that 
{Y(t) : t > 0} is a Markov process with state space F , evolving according to a transition 
semigroup (Q t ; t > 0) and with initial distribution v . Suppose further that L is a Markov 
transition kernel from E to F, such that [iL = v and the intertwining PtL = LQ t holds. 
Now let f : E — > G and g : F — > G be maps into a third state space G, and suppose that 

L(x, •) is carried by {y G F : g{y) = f(x)} for each x G E. 

Then we have 

{f(X(t)) :t>0}± {g{Y{t)) : t > 0}, 
in the sense of finite dimensional distributions. 

Proof of Lemma^ For any bounded function a on G let Tia be the function ao f defined 
on E and let I^a be the function a o g defined on F. Then it follows from the condition 
that L(x, ■) is carried by {y G F : g(y) = f(x)} that whenever h is a bounded function 
defined on F then 

L(T 2 a x h) = r x a x Lh, (2.1) 

which is shorthand for J L(x , dy)T 2a(y)h(y) = Via x Lh. For any bounded test functions 
ao,ai, • ■ ■ , a n defined on G, and times < t\ < ■ ■ ■ < t n , we have, using the previous 
equation and the intertwining relation repeatedly, 

E[a o (^(y(0)))a 1 (^(y(t 1 ))) . . . a n {g(Y{t n )))\ 
= v(T 2 a x Q tl (T 2 ai x Q t2 ~ tl (- ■ ■ (r 2 a„-i x Qt„- 4n _ i r 2 a n ) •••))) 
= /iL(r 2 a x Q tl (r 2 «i x Q t2 - tl (• • • (r 2 a„_i x (5 tn _ tn _ i r 2 a; n ) •••))) 
= /i(ria x P^iYxai x P t2 - tl {- ■ • (r X Q; n _i x P tn _i n _ i ria„) •■■))) 
= E[a (/(X(0)))a 1 (/(X(t 1 ))) . . . a n (f(X(t n )))} (2.2) 

which proves the equality in law. □ 

We let (Y(t) : t > 0) be the process Y of n reflected Brownian motions with a wall 
introduced in the previous section. It is clear from the construction (11.101) that the pro- 
cess Y is a time homogeneous Markov process. We denote its transition semigroup by 
(Qt]t > 0). It turns out that there is an explicit formula for the corresponding den- 
sities. Recall (f) t (z) = y 2 =7 e ~ z2 ^ 2t ' ) - Let us define 4>t\y) = ^k4>t(y) for k > and 

<pt k \y) = (-i) fc XT bj$^M*)dz for k> 1. 
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Figure 1: The set IK. The triangle represents the intertwining relations of the variables z 
and the vertical line on the left indicates zf k+1 > 0, see (12. 5ft . (12. 6ft . The set of variables on 
the bottom line is denoted by b(z) and the one on the upper right line by e(z). 

Proposition 5. The transition densities qt{y,y') from y = (yi, . . . ,y n ) at t = to y' = 
(y[, . . . , y' n ) at t of the Y process can be written as 

Qt(y, y') = det{aij(yi, yj)}i<<j<n (2.3) 

where ajj is given by 

auto, y') = (-irV?~% + y) + (-iy + ^t l \y - y')- (2-4) 

The same type of formula was first obtained for the totally asymmetric simple exclusion 
process by Schiitz [13]. The formula for the Z process was given as a Proposition 8 in [15], 
see also [12]. 

Proof of Proposition^ For a fixed y', define G t (y,t) to be (12.31) as a function of y and t. 
We check that G satisfies (i) the heat equation, (ii) the boundary conditions J^| yi=0 = 0, 
f^|yi=2/i-i = °> i = 2, 3, . . . , n and (iii) the initial conditions G(y, t = 0) = YYi=i ~~ V'd- 
(i) holds since 4>[ k \y) for each k satisfies the heat equation, (ii) follows from the 

relations, ^a l:j (y,y% = o = <j>t\y') + (-l) i+ Vt (~y') = and ^ a ij{y^') = -a*-u(l/,y')- 
For (iii) we notice that the first term in (12 .4p goes to zero as t — > for y, y 1 > and the 
statement for the remaining part is shown in Lemma 7 in [15]. □ 

For n = 2m, resp. n = 2m — 1 we take (X(t),t > 0) to be Dyson Brownian motion of 
type C, resp. of type D. The transition semigroup (Pt, t > 0) of this process is given by 

(O. 

Let IK denote the set with n layers z = (z 1 , z 2 , . . . , z n ) where z 2k = [z\ k , z 2k , . . . , zf fc ) e 
R^, z 2k ~ l = (z 2fc ~\ zf~ x , . . . , z 2 ^ v ) E and the intertwining relations, 

z 2k-l < z 2k < ^fc-l < z 2k < < ^k-l < z 2fc (2 5) 

and 

< ^ 2fc+1 < ^ 2fc < ^ 2 2fc+1 < 4 fc < . . . < ^ 2fc < ( 2 - 6 ) 
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hold (Fig. [T]). Let n = 2m or n = 2m — 1 for some integer m. We define a kernel L° 
from E = {0 < x\ < . . . < x m } to F = {0 < y\ < . . . < y n }. For z G K, define 
b( z ) = z 11 = (z?, . . . , z£) G E, e(z) = {z\, zl zl 4, . . . , z£) E F and K(x) = {z E 
K; b(z) = x G E}, K[y] = {zGK; e(z) =y G F}. The kernel L° is defined by 

L^(x)= / L°(x,dy)g(y)= [ g(e{z))dz. (2.7) 
Jf Jk(x) 

where the integrals are taken with respect to Lebesgue measure but integrations with 
respect to z on the RHS is for b(z) = x fixed. 

The function h defined at (11.61) is equal to the Euclidean volume of K(x). Consequently 
we may define L to be the Markov kernel L(x, dy) = L°(x, dy)/h(x). In the remaining part 
of this section we show 

Proposition 6. 

LQ t = P t L. (2.8) 

Now if we apply Lemma H] with f(x) = x m , g(y) = y n and the initial conditions starting 
from the origin we obtain fl!.13j) . 

Proof of Proposition^ The kernels Pt{x, •) and L(x, •) are continuous in x. Thus we may 
consider x in the interior of E, and it is enough to prove 

(L°Q t )(x,dy) = (P t L°)(x,dy). (2.9) 

From the definition of the kernel L°, this is equivalent to showing 

/ q t (e(z),y)dz= [ p° t (x,b{z))dz (2.10) 

JK(x) JK[y] 

where q t and p° are densities corresponding to Q t and P t °. Integrations with respect to z 
are on the LHS with b(z) = x fixed and on the RHS with e(z) = y fixed. 

Let us consider the case where n = 2m. Using the determinantal expressions for q t 
and p° we show that both sides of (12.101) are equal to the determinant of size 2m whose 
matrix element is a 2 ij(0,yj) for 1 < % < m, 1 < j < 2m and a 2m ,j(%i-m,yj) f° r 
m + 1 < % < 2m, 1 < j < 2m. 

The integrand of the LHS of (l2~T0l is 

q t (e(z),y) = det{a i) j(e(2;) i , yj)}i<i,j<2m ( 2 -H) 

with b(z) = x. We perform the integral with respect to z 1 , . . . , z 2m ~ x in this order. After 
the integral up to z 21-1 , 1 < / < m, we get the determinant of size 2m whose matrix 
element is a 2 ij(0,yj) for 1 < i < I, a2i^{zf_ l ,yf) for / + 1 < i < 21 and a it j(e(z)i,yj) for 
21 + 1 < i < 2m. Here we use a property of a^j, 

POO 

a id (y,y')= a,i-i d (u,y')du, (2.12) 
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and do some row operations in the determinant. The case for I = m gives the desired 
expression. 

The integrand of the RHS of fl2~T0D is 



Pt(x,z 2m ) = det(a 2m ,2m(a;i,2;- m ))i<ij<r 



(2.13) 



with the condition e(z) 



We perform the integrals with respect to 



y 2m 
'1 i 



J2m 



2m- 1 



%m-\ )>•••) z ii z i m this order. We use properties of 



■3i 



%j(y,y') 



a itj+1 (y,u)du, 



a 2 i,2j(x,0) = 0, a 2 i,2i-i(0,y) = 1, a 2i j(0,y) = 0, 2i < j. 



(2.14) 
(2.15) 



After each integration corresponding to a layer of K we simplify the determinant using 
column operations. We also expand the size of the determinant after an integration corre- 
sponding to (zf , . . . , zf_i) for 1 < I < m, by adding a new first row 

(M__ 1 i,o 1 o__ 1 o) = 

l 2m-2Z+l 

(oa,2i-i(0, Z?' 1 ), . . . , aa,2i-i(0, zf 1 ' 1 ), a 2lt2l (Q, e(z) 2l ), . . . , a 2 i ;2m (0, e{z) 2m ))) (2.16) 

together with a new column. After the integrals up to (z\ _1 , . . . , zfi[ l ) have been per- 
formed, we obtain the determinant of size 2m — I + 1, 



O2(/+j-l),2(/-l)(0, ^ ^) «2(/+i-l),j+Z-l(0, e(z)j + i-i) 



Ct>2m,2(l-l)(%i-m+l-li z j ) a 2m,j+Z-l (^i-m+Z-1 , e (^)j+i-l) 



(2.17) 



Here l<2<m — Z + l (resp. m — l + 2<i< 2m — I + 1) for the upper expression (resp. 
the lower expression) and 1 < j '• < I — 1 (resp. I < j < 2m — l + l) for the left (resp. right) 
expression. For I = 1 this reduces to the same determinant as for the LHS. 

The case n = 2m — 1 is almost identical. Similar arguments show that both sides of 
( 12. 9p are equal to a determinant size 2m — 1 whose matrix element is a 2i j(0,yj) for 
1 < % < m — 1, 1 < j < 2m — 1 and a,2m-i,j( x i-m+ii Vj) f° r m + 1 < i < 2m — 1, 1 < j < 
2m -1. □ 



3 Proof of proposition [3] 

Using (11.101) repeatedly, one has 



Y n (t) 



sup 

0<ti<...<t n <t 



i+lj 



Bi(U)) 



(3.1) 



i=i 
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with t n+ i = t. By renaming t — t n _ i+ i by t{ and changing the order of the summation, we 
have 

n 

Y n (t)= sup J2{B n _ i+1 (t-t i+1 )-B n .. i+1 (t-t i )). (3.2) 

0<ti<...<tn<t ~q 

Since B^s) := B n _ i+X {t) - B n _ i+1 (t - s) = B^s), 

n 

Y n (t) = sup YYfliCtj) - Bi(t - U-i)) = sup Z n (t). (3.3) 

0<h<...<t n <t 0<s<t 
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